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Abstract 

Markovian diffusion processes yield a system of conservation laws which cou- 
ple various conditional expectation values (local moments). Solutions of that 
closed system of deterministic partial differential equations stand for a regular 
alternative to erratic (irregular) sample paths that are associated with weak solu- 
tions of the primordial stochastic differential equations. We investigate an issue 
of local characteristics of motion in the non-Gaussian context, when moments of 
the probability measure may not exist. A particular emphasis is put on jump- 
type stochastic processes with the Ornstein-Uhlenbeck-Cauchy process as a fully 
computable exemplary case. 

1 Local characteristics of the Brownian motion 

Lets us consider a Markov process Xt in . We can fully characterize an associated 
random dynamics by means of a transition density p(?/, s,x,t) with < s < t and an 
initial density po{x) = p{x, to),0 < to. For the stochastic process to be properly defined 
we impose the so-called stochastic continuity condition: 

limtis-— I p{y,s,x,t)dx = (1) 

t — S J\y—x\>e 

to be valid for arbitrary (every) e > and for almost every y & R^. That is known to 
imply that along a sample path there holds /im(jsP[u; : \Xt{uj)—Xs{uj)\ > e] = 0. This 
condition needs be respected by both diffusion and jump-type propagation scenarios 
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and is known to be necessary but still insufficient for the Markov process to have 
continuous trajectories. 

Various characteristics of random motion can be directly associated with proper- 
ties of sample trajectories (consider for concretness the iterated logarithm laws in the 
Brownian motion). Those characteristics of random dynamics, that are free of irregu- 
larities of sample paths, need to involve some forms of averaging. That refers either to 
global averaging present while evaluating moments of a probability measure (provided 
they exist) or to the evaluation of local moments of that measure. 

Given a transition probability of a Markov process, we can attempt to derive a 
number of local expectation values, like e.g. the forward drift of the process: 

b{x, s) = limtis-^ I {y- x)p{x, s, y, t)dy (2) 

t — S J\y—x\<& 

and the diffusion (coefficient if a constant) function 

D{x,s) = limtis-^ [ {y - xfp{x,s,y,t)dy (3) 

t — S J\y—x\<S 

where the S cutoff is needed to guarantee a convergence of the integral. In principle, 
the drift and diffusion functions should take values independent of the 5-cutoff and/or 
exist when the cutoff is removed (once we let S go to oo). However, for the jump- 
type processes, local moments (e.g. drifts and diffusion functions) in the absence of a 
cutoff are generically nonexistent, while in a cutoff version they display an explicit and 
nontrivial 5-dependence. 

As yet, our discussion extends to both continuous and discontinuous processes. Let 
us explicitly reveal when a diffusion process enters the game. By assuming that there 
holds: 

limtisz / \y - x\'^+'^p{x,s,y,t)dy ^0 (4) 

t — s Jr 

for any 7 > 0, we set a sufficient condition for the continuity of sample paths (to hold 
true almost surely) , that also allows to remove the previous 5-cutoff from local momet 
formulas. The resulting process is a Markovian diffusion process with the forward drift 
b{x,s) and diffusion function D{x,s). Our futher discussion will be carried under a 
simplifying assumption that a diffusion function D{x, s) actually is a constant which 
we denote D. 
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Our previous considerations, when specialized to Markovian diffusion-type pro- 
cesses, can be cast in another form. Namely, we can depart from a formal infinitesimal 
version of a stochastic differential equation for a random variable Xf — X{t) taking 
values in R^: 

dXt = b{X(t),t)dt + V2DdWt (5) 

where Wt = W{t) is a Wiener process and b{x, t) stands for the forward drift of the 
diffusion-type process Xt. 

If we assign a probability density po{x) with which the initial data Xo — X{0) for the 
stochastic differential equation are distributed (weak solutions enter the scene), then 
the emergent Fick law would reveal a statistical tendency of particles to flow away 
from higher probability residence areas. This feature is encoded in the corresponding 
Fokker-Planck equation (equivalently, a continuity equation): 

dtp = -V • (vp) (6) 

where a diffusion current velocity is 

(7) 

Clearly, the local diffusion current (a local flow that might be experimentally ob- 
served for a cloud of suspended particles in a liquid) j = vp is nonzero in the nonequi- 
librium situation and quantifles a non-negligible matter transport which occurs as a 
consequence of the Brownian motion, on the ensemble average. 

It is interesting to notice that the local velocity fleld v{x, t) obeys the natural (local) 
conservation law, which we quite intentionally pattern after the moment identities 
(hierarchy of conservation laws) valid for the Boltzmann and Kramers equations. The 
pertinent momentum conservation law directly originates from the rules of the Ito 
calculus for Markovian diffusion processes, and from the flrst moment equation in the 
diffusion approximation of the Kramers theory: 

dtV + {vV)v^Vin-Q) . (8) 

The general form of an auxiliary potential fl{x,t) reads: 

n(x,t)=2D[dt<j> + l{^ + V-b)] . (9) 
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and can be interpreted as a constraint equation for the admissible forward drift b{x, t) — 
2D'V(f){x, t). It is useful to mention that the above momentum conservation law may be 
regarded to account for effects of external (conservative) volume forces — V(— Q) which 
are determined by a priori arbitrary, but bounded from below, continuous function 
Q{x^t). Then, the allowed forward drift must be disentangled from the Ricatti-type 
equation (9) as a function of ft. 

In the above there appears a contribution from a probability density p-dependent 
potential Q{x,t). It is is given in terms of the so-called osmotic velocity field u{x,t): 

Q{x,t) ^^u^ + -u (10) 

u{x,t) = DV lnp{x,t) 

and is generic to a local momentum conservation law respected by isothermal Marko- 
vian diffusion processes. Notice that in case of the free Brownian motion (admitted, if 
we set = 0), we would have v{x,i) = —u{x,t) for all times. An equivalent form of 
the pressure- type potential Q is Q — 2D — 175" • 

It is interesting to observe that a gradient field ansatz for the diffusion current 
velocity {v — VS): 

dtp = -V • [{WS)p] (11) 

allows to transform the momentum conservation law of a Markovian diffusion process 
to the universal Hamilton-Jacobi form: 

n^dtS + -\vs\^ + Q (12) 

where Q{x,t) was defined before. By applying the gradient operation we immediately 
recover the previous local momentum conservation law. 

Notice that the above Hamilton-Jacobi- type equation is sensitive to any additive 
(constant or time-dependent) modification of the potential Q. In the above, the con- 
tribution due to Q is a direct consequence of an initial probability measure choice for 
the diffusion process, while Q alone does account for an appropriate forward drift of 
the process. 

The simplest realisation of the outlined theoretical framework can be provided 
by invoking a standard Ornstein-Uhlenbeck process. Namely, let us consider an Ito 
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equation (in its symbolic differential version) for infinitesimal increments of the velocity 
random variable, exhibiting the systematic frictional resistance: 

dV{t) = -pV{t)dt + pV2DdW{t) (13) 

where W{t) is the normalized Wiener process. One can easily infer the corresponding 
second Kolmogorov (Fokker- Planck) equation 

dtp{vo, V, t) = f]DA^ + pV^ ■ [vp{vo, V, t)] (14) 

for the transition probability density of the time homogeneous process in the velocity 
space alone. The pertinent transition probability density reads: 

Let us consider an instaiitaucous velocity Vt = that has been achieved in the 
course of the Ornstein-Uhlenbeck random evolution beginning from a certain Vq = "^o- 
We can evaluate a conditional expectation value (local mean with respect to the law 
of random displacements) over all randomly accessible velocities V{t + At) = at a 
time t + At, At > 0. That determines the forward drift of the process: 

b{v, t) — liiriAtioi j v' p{v, v', At) dv' — v] — —(3v (16) 

and thus provides us with an information about the mean tendency of the dynamics on 
small (but not too small if compared to the relaxation time (i^^) scales. Analogously, 
we can derive a diffusion function for the Ornstein-Uhlenbeck process which is constant 
and equals 20(3"^. 

To arrive at local conservation laws, an additional input of an initial probability 
density po{v) of Vq is necessary. To that end, one may choose an asymptotic stationary 
(invariant, Maxwell-Boltzmann) density of the Ornstein-Uhlenbeck process: po{v) — 

( _i]2_)l/2 trmPl 

2 Prom Poisson probability law to jump - type pro- 
cesses 

A random variable X taking discrete values 0,y,2y,3y, with y > is said to have 
Poisson distribution V{X), A > with jump size y, if the probability oi X = ky is given 
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by P{X — ky) — ^exp{—X). The characteristic function of V{\) reads: 

\k 

E[exp{ipX)\ = exp[\{e'Py - 1)] = E^e"^ e'^^^^ = E^e^'^^f^) P{X = ky) (17) 

K ! 

and the first moment of the probabihty measure equals E[X] = A. Notice that 
P{X = 0) = exp{—X), hence the numerical value of A > fixes the no-jump proba- 
bility. For the Poisson random variable with values b + ky,k = 0, 1, ...we would get 
E[exp{ipX)] = exp[ibp + A(e*^^ — 1)]. 

If we consider n independent random variables Xj,l < j < n such that Xj has 
Poisson distribution V{\j) with jump size yj, then a new variable X can be introduced 
by means of the distribution of Xi + ... + X„ whose characteristic function reads 

E[exp{ipX)] = exp[TP^^^\{e'Py^ - 1)]. (18) 

The exponent in the above might include an additional term ipE^bj if nonrandom 
shifts of each jump kyj by bj were allowed. 

We can admit not only jumps of fixed magnitudes yi, ■■■,yn but also jumps covering 
an arbitrary range in R+. Let the distribution function of the magnitude of the jump 
be P{x < y) — iJi{y). In this case we set 

E[exp{ipX)] = exp[ f {e'^^ - l)ii{dy)] (19) 

assuming that the integral in the exponent exists. (Notice that the previous formula 
is recovered, if we choose ii{dy) = Ti^^iXj6{y — yj) dy.) 

For any Borel set A <Z R bounded away from the origin, the random variable X^ 
representing jumps bounded by A, gives rise to a characteristic exponent /^(e*™ — 
l)dfi{y), and the expected number of jumps of size bounded by A is equal to 

IJ,{A). We can interpret that in terms of jumps of different sizes that are mutually 
independent. Jumps whose size is bounded by [y,y + Ay), Ay <^ 1, do contribute a 
Poisson component with exponent function approximately equal to (e*^^ — l)ii{[y, y + 
Ay)). 

Let us consider an expression for a characteristic function of a probability measure of 
a certain random variable X that is given in the form E[exp{ipX)] = exp[—F{p)] where 
p e and for — oo < p < +00, F = F{p) is a real valued, bounded from below, locally 
integrable function. If F{p) satisfies the celebrated Levy-Khintchine formula, then the 
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pertinent measure is positive and we may introduce positivity preserving semigroups, 
together with the induced (Markovian) stochastic processes. Let us concentrate our 
attention on the integral part of the Levy-Khintchine formula, which is responsible for 
arbitrary stochastic jump features. In that case, F{p) takes the form: 

F{p) = - r^[exp{ipy) - 1 - -^]„{dy) (20) 

where i^idy) stands for the so-called Levy measure on R^. 

A generic feature of jump-type processes is that they admit jumps of arbitrarily 
small size (without any lower bound) and some care is necessary when evaluating 
contributions from a close neighbourhood of the origin. One obvious way to bypass this 
problem amounts to an e-cutoff which allows to neglect "too small" jumps. Namely, for 
a probability law with the characteristic exponent — F(p), we can consider its restriction 
to upward jumps of size exceeding a given lower bound like e.g. to all y > e > 0: 

4>tip) = I - 1 - T^lK^y) = / V"^ - ^My{y) - ipK (21) 

Clearly, we deal here with a random variable of the type considered before, and we 
can try to isolate contributions from jumps of the size [y, y^ Ay) by coarse-graining a 
Borel set A of interest. A formal exploitation of 

u{dy) = J:]=,Xj5{y - yj) dy (22) 

gives rise to 

E[exp{ipX)] = exp[E]^^ [Xj {e'^^^ - 1) - ip r^^] . (23) 

Further specializing the problem we shall consider Levy measures that obey the 
spatial reflection property u{—dy) = —v{dy). Then, we can readily extend our dis- 
cussion to jumps of all sorts in R}, i.e. y can take values in both i?+ and with 
the only restriction to be observed that \y\ > e > 0. Notice that we shall deal with 
two distinct types of jumps, either positive or negative, with no common jump point 
for them. This fact means that they are independent components of the more general 
random variable: 

Up) -I [e'^-lHdy)-ipK (24) 

J\y\>e 
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where (choose ^{dy) — for concretness) the deterministic term identically vanishes 
in view of 

6, = ht + b:= [ Y^A^y) + / vh^^'^'^y^ ^ o- ^^s) 

Jy>e l + y^ Jy<-e 1 + 

In the previous steps we have indirectly exploited a defining property of infinitely di- 
visible probability laws: if exp(/){p) is a characteristic function of a given probability dis- 
tribution, then [exp4){p)Y = exp[t(f){p)], t > is likewise a characteristic function of an 
infinitely divisible probability law again. This feature extends our discussion to stochas- 
tic jump and jump-type processes (time homogeneous with independent increments). 
Obviously, for such processes E[exp{ipX (t))] = exp[t(l){p)] while EA[X{t)] = tu{A), 
and our previous arguments retain their validity with respect to 

E[exp{ipX{t))]e^exp[t(f)e{p)]^exp[t [ {e'^ - l)u{dy)]. (26) 

Coming back to the Lcvy-Khintchine formula, let us replace a function F{p) by an 
operator acting in a suitable domain according to the recipe: F[p) H = F[p) where 
p = — zV. We easily learn that for times t > there holds 

[exp{-tH)]f{x) = [exp{-tF{p))f{p)]'^{x) (27) 

where the superscript V denotes the inverse Fourier transform and / stands for the 
Fourier transform of a function /. 

If we set Pt = -^[exp{—tF{p)Y, then the action of exp{—tH) can be given in terms 
of a convolution: exp{—tH)f — f * Pt, where (/ * g){x) :— Jugix — z)f{z)dz. Clearly, 
there holds: 

dtp{x,t)^-{Hp){x,t)^ (28) 
dpe{A, t)= dx[ / [xa{x + y)- XA{x)]u{dy)]pe{x, t) 

JR J\y\>e 

which displays a generic Master-equation form. Indeed, we have here dtPe{A,t) = 
Jj^q^{x,t, A)p^{x,t)dx where is interpreted as the jump intensity. 

Let us however emphasize that the above simplification occurs only in the \y\ > e > 
jumping size regime. The real role of two integral terms in the expression for is 
to compensate the divergent contributions from the Levy measure when the principal 
value integral e — > limit is considered; then the standard jump process theory does 
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not literally apply since arbitrarily small jumps are admitted. Anyway, those two terms 
are irrelevant if we assume an e > cutoff, irrespective of how small the chosen (and 
fixed) e is. 

The best known example of the stable probability law that is compatible with the 
above definitions is provided by the classic Cauchy density which will be our reference 
probability law in below. Let us focus our attention on that case which is specified by 
F{p) — \p\. The corresponding semigroup generator H — \V\ is a pseudodifferential 
operator. The associated kernel pt is a transition density of the jump-type Cauchy 
process, which is a solution of a pseudodifferential Fokker-Planck equation: 

dtp{x,t) = -\V\p{x,t). (29) 

The pertinent probability density reads p{x,t) = ^-fj^ and the corresponding 
space-time homogeneous transition density (e.g. the semigroup kernel function) is: 

'^^l ^ ^ p^y^ ^' ^ l it-srllx-yr ^ 

0< s<t 

{exp[ipX(t)]) :— / exp{ipx)p{x,t) dx — exp{—\p\t). 
Jr 

The characteristic function of p{y,s,x,t) for y,s fixed, reads exp[ipy — \p\{t — 
s)], and the Levy measure needed to evaluate the Levy-Khintchine integral reads: 
iy{dy) ■.= Hmtio[^^piO,0,y,t)]dy=^. 

The Cauchy process belongs to the category of jump-type processes, where apart 
from the long jumps-tail (no fixed bound can be imposed on their length) which is 
the reason of the nonexistence of moments of the probability measure, sample paths of 
the process may have an infinite number of jumps of arbitrarily small size. By general 
arguments, pertaining to the space De[0,oo) of right continuous functions with left 
limits (cadlag), both in the finite and infinite time interval the number of jumps is at 
most countable. 

An introduction of the e-cutoff which eliminates small jumps is fed up by the physi- 
cal intuition. An approximation of the jump-type process should in principle be possible 
in terms of more traditional jump processes which involve a finite number of jumps in 
a finite time interval. This radical approximation assumption is usually achieved by 
giving a characterisation of the affiliated Markovian jump-type processes in terms of 
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approximating (convergent) families of so-called step processes. The step processes are 
not yet the jump processes of the standard daily experimental evidence. They need to 
have no accumulation points of jumps in a finite time interval: in that case the number 
of jumps is finite on each finite time interval and between jumps the sample path is 
constant. 

Let us recall that the operator |V| acts as follows: 



By turning back to the pscudodifferential Fokker-Planck equation with the e-cutoff 
implemented, let us introduce an operator IVI^: 



It suffices to replace / by to arrive at the right-hand-side of the previously defined 
Fokker-Planck equation for the step process approximant of the Cauchy process. Its 
generator is just |V|e. 

While mitigating the " arbitrarily small jumps" problem the e-cutoff does not remove 
all obstacles related to the Cauchy process. Indeed, at the first glance the situation 
looks deceivingly simple, because on a finite time interval there can be at most finitely 
many points t e [0, T] at which the jump size exceeds a given positive number. In 
view of that, supt^Yi,T\ \^t \ < °° where XI stands for the e-bounded Cauchy process 
(the same argument extends to the unrestricted Cauchy process Xf). However, there 
must be no fixed upper bound for the size of jumps (except for being finite), since a 
stochastically continuous process with independent increments having, with probability 
1, no jumps of size exceeding a certain constant 5, would possess all moments. That is 
certainly not the case for the Cauchy process, which is known not to have any moments. 

Hence imposing or not imposing an upper bound on the jump size (call it a 5 
cutoff) is another critical issue that hampers a reliable approximation of the jump - 
type process in terms of experimentally verifiable jump processes whose jump size is 
bounded both from below and from above. That derives in part from to the resolution 
limitations of realistic experimental arrangements (any experimental data collection 
and any computer simulation/experimentation have built-in lower and upper jump 
size bounds), and in part from the fact that all observations are carried in finite time 



V|/(x) = -- / [f{x + y)-f{x) 



yVf{x) ^dy 
y2- 



(31) 




(32) 
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on systems of finite spatial extension. Let us point out that by imposing both e and 
5 cutoffs on the Cauchy process, we would reduce the problem to the standard jump 
process. 

Obstacles related to heavy tails of the probability distribution can be visualized by 
making computer experiments for the converging (in fact, diverging) variance test. 
Namely, having given a sample of jumps determined by the Cauchy distribution: 
Xi, X2, ...Xn, with 1 < i < n,one can form a statistics based on first i "observa- 
tions" and ask for the behaviour of the averages with respect to i. First we need the 
i-th mean value: X — j J2k=i ^^xt the i-th variance: Sf — Sfc=i(^fe — Xi^. 

Plots of Xi and Sf against i would show up a fairly irregular behaviour and definite 
non-convergence signatures for large values oi i < n. 

3 Local characteristics of the Ornstein-Uhlenbeck- 
Cauchy process 

Although we directly refer to the specific Cauchy stochastic process, in fact we stay in 
a much broader setting of so-called Levy flights. Generically, the variance and higher 
cumulants of those processes are infinite (nonexistent). There is also physically more 
singular subclass of such processes for which even the first moment (mean value) is 
nonexistent. That is true for the Cauchy process. Thus we need to relax the limitations 
of the standard Gaussian paradigm: we face here a fundamental problem of establishing 
other means (than variances and mean values) to characterize statistical properties of 
Levy processes. 

Specifically, if a habitual statistical analysis is performed on any experimentally 
available set of frequency data, there is no obvious method to extract a reliable infor- 
mation about tendencies (local mean values) of the random dynamics. Nonexistence 
of mean values and higher moments may also be interpreted as the nonexistence of 
observable (e.g. mean, like drifts or local currents) regularities of the dynamics. More- 
over, as we have learned before, the jump-type processes usually admit arbitrarily 
small jumps (with no lower bound) and finite, but arbitrarily large jump sizes (with 
no upper bound). Any laboratory experiment or computer simulation would involve 
both the lower (coarse-graining) and upper bound on the jump size. Mathematically, 
that sets (as suitable) the framework of standard jump processes for which the central 
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limit theorem is known to hold true in its Gaussian version (even if we account for a 
possible abnormally slow convergence to a Gaussian). Therefore, there is no clear-cut 
procedures allowing to attribute an unambiguous statistical interpretation in terms of 
Levy processes to given phenomenological data. Moreover, no realistic formulation of 
a fluctuation-dissipation theorem is possible in that case (nonexistence of variances) 
which pushes us away from any conceivable thermal equilibrium framework. 

The starting point for Ornstein and Uhlenbeck (cf. the previous Section) was the 
dissipative Langevin equation 

^ = -XV{t) + A{t) (33) 

where V{t) is a random variable describing the velocity of a particle, A > is a 
dissipation constant, and A(t) is another random variable whose probabilistic features 
are determined by the probability distribution of V{t), which is assumed to satisfy a 
concrete law when t ^ oo. Because V{t) may have no time derivative, the Langevin 
equation was soon replaced by another one, the stochastic differential equation, namely 

dV(t) = -XV(t)dt + dB{t), y(0) = vo (34) 

which received a rigorous interpretation within the framework of stochastic analysis. 
In the case when the probability distribution of V{t), t oo, is the Maxwell one, B{t) 
must be a Gaussian (in fact Wiener, Bt — Wt) process, and then we end up with a 
classical Ornstein-Uhlenbeck process. 

Now, we shall discuss properties of the process V^(t), in the case when B = {B{t))t>Q 
is the Cauchy process instead of the traditional Wiener one. By straightforward inte- 
gration we obtain that for i > s 

t 

V{t) = e-^(*-^V(s) + e-^' J e^^dB{r). (35) 

s 

There are a number of (eqiiivalent) procedures to deduce a probability density 
of the process V{t) from the Cauchy increments statistics. One may follow a direct 
probabilistic route which, upon assuming that a characteristic function of the Cauchy 
probability measure reads 
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with the choice of ip{p) — (T^|p| (before, we have used cr^ = 1), leads to the transition 
density: 

pUu,v) = piv(t) = urns) = .] = + 

where a^{t - s) = ^(1 - e-^(*-^)). 

Since V{0) = Vq, the probabihty density of V{t) is given by P[V{t) = v] — 

1 T^t) 



TT (i;-'i;oe-^*)2 + <T4(t)- 

Now, we shall demonstrate an important property (mentioned before in connection 
with the Ornstein-Uhlenbeck process) of the so-called stochastic continuity which is 
a necessary condition to give a stochastic process an unambiguous status. Namely, 
we need to show that for any e > the following equation is satisfied limF[|y(t) — 



> e] = 0. This equation is equivalent to lim / pt{u,v)du = 0. Because of 

*^0|«-D|>e 
;(1 - e~^') I , e-v(l-e 

"5^(t) ^ arctan 



/ pt{u\v)du = 1 — ;^ [arctan — - + arctan^ — ^%(t) — ''\ and remembering 



\u—v\>e 

that cr'^{t) — y-(l — e"'^*), the stochastic continuity property does follow. 

The nonexistence of moments of the probability measure in case of the Cauchy 
process leads to straightforward difficulties, since the standard local characteristics of 
the diffusion-type process like the drift and the diffusion function (or coefficient) seem 
to be excluded in the present case. However, for the considered Ornstein-Uhlenbeck- 
Cauchy process, the notion of the forward drift of the process proves to make sense. 

Since we know the Markov transition function pt_s(M, v), t > s, for the process Vt, 
we can exploit our experience with diffusion processes and say that the process Vf has 
a drift (in fact, forward drift) if the following limit 

lim / (u — v)pt-s{u,v)du (38) 

t — s J 

\u—v\<5 

does not depend on the choice of 5 > 0. If so, then its value depending only on (v, s) 
we denote by b{v, s) and call it the drift coefficient. Clearly, if p is homogeneous in 
time, then the drift coefficient depends only on the variable v. Let us emphasize that 
in the above definition we do not require the process Vt to have finite moments. 
We claim that the jump-type Markov process V{t) has a (forward) drift which reads 
b{v) = —Xv. Indeed, by first evaluating the indefinite integral 

1 f. a'^(t)du 
I ^ - (u 
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and substituting z — u — ve we get 

r zdz V _ f a\t)dz 

n J + a\t) TT^ ^Jz^ + a^{t) 

- '^^^^hogiz'' + a\t)) + -(e-^* - l)arctan(- ^ 



Hence 



/ ^ ^\og[{u - ve-^'f + a\t)] + -(e-^* - 1) arctan['' '''' 



27r ' ' tt' ' ' (72(i) 

and consequently the limit 

t^o t 



lim - /!"=!+! 



linil^(log[(^; + e - ^;e-^*)2 + (7^(i)] - log[(t; - e - ve-^'f + 

Iv , _xf ;V + e — ve~^* , .V — e — ve"'^* , 
+ lim — e — 1 arctan -— — arctan r— 

'V .71 TT. 

= - A-(- + -) = -A^ 

TT Z Z 

exists and is e-independent. This is the forward drift of the process V{t) which proves a 
consistency of the derived transition probability density with the stochastic differential 
equation for the process Vt — V{t) . 

It is well known that for Markovian diffusion processes all local characteristics of 
motion (conditional expectation values that yield drifts and variances) are derivable 
from transition probability densities, supplemented (if needed) by the density of the 
process. We have demonstrated that, in the non-Gaussian context, the nonexistence 
of moments does not necessarily imply the nonexistence of local characteristics (drifts) 
of the process. 

However, the situation becomes uncomfortable once we attempt to evaluate another 
local moment. Namely, in the present case there holds: 

1 r 2(j^ 
lirritis / (u — v)'^pt(u,v)du — S (39) 

t — S J\v-u\<S TT 

i. e. an explicit cutoff d (upper bound on the size of jumps) persists in this formula 
and there is no way to remove that jump size restriction from the formalism, unless we 
wish to get the divergent integral. 
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This property is a clear indication that a convergence to a Gaussian might always 
be expected if the Ornstein-Uhlenbeck-Cauchy process is approximated (we disregard 
an issue of how good that approximation is) by means of jump processes with an upper 
and lower bound on the jump size. In that case both the mean and variance would 
exist for the approximating process. In particular, the central limit theorem would 
work as usual for the (e, 5)-jump process approximation of the Cauchy process Bf. 
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